The Astrophysical Journal, 683:L207-L210, 2008 August 20 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



oo 
O 
O 
(N 

D 
(N 

O 



> 

O 
oo 
O 



X 



EVIDENCE OF SHOCK-DRIVEN TURBULENCE IN THE SOLAR CHROMOSPHERE 

K.P. Reardon^-^, F. Lepreti^, V. Carbone'*"', and A. Vecchio^"' 

Received: 2008 March 18; accepted 2008 July 11; 2008 August 12 

ABSTRACT 

We study the acoustic properties of the solar chromosphere in the high-frequency regime using a 
time sequence of velocity measurements in the chromospheric Call 854.2 nm line taken with the 
Interferometric Bidimensional Spectrometer (IBIS). We concentrate on quiet-Sun behavior, apply 
Fourier analysis, and characterize the observations in terms of the probability density functions (PDFs) 
of velocity increments. We confirm the presence of significant oscillatory fluctuation power above the 
cutoff frequency and find that it obeys a power-law distribution with frequency up to our 25 mHz 
Nyquist limit. The chromospheric PDFs are non-Gaussian and asymmetric and they differ among 
network, fibril, and internetwork regions. This suggests that the chromospheric high-frequency power 
is not simply the result of short-period waves propagating upward from the photosphere but rather is 
the signature of turbulence generated within the chromosphere from shock oscillations near the cutoff 
frequency. The presence of this pervasive and broad spectrum of motions in the chromosphere is likely 
to have implications for the excitation of coronal loop oscillations. 
Subject headings: Sun: chromosphere, turbulence, shock waves. Sun: photosphere 



1. INTRODUCTION 

The role played by high-frequency acoustic oscillations 
in the heating of the solar chromosphere and corona has 
been addressed in numerous studies. Because acoustic 
oscillations above the cutoff frequency freely propagate 
upward, they have long been considered candidates for 
transporting to the outer solar atmosphere some of the 
abundant mechanical energy generated below the ph oto- 
sphere (see review bv lNarain fc Ulmschneidalll996l ). 

The primary test for this theory has been to measure 
the high-frequency oscillations in the photosphere and to 
estimate whether they contain sufficient acoustic flux to 
balanc e the chromospheri c and coro nal losses. Fo r ex- 
ample, ImnMnbirieE^ (|200l) andlAndi^ (|2007f) have 
looked for power in velocities and in tensities measured in 
photospheric lines. Other authors f de Wiin et al.l [20051 : 
iDeForest 2004 ) have used Transition Region and Coro- 
nal Explorer [TRACE) UV continuum images to mea- 
sure power spectra of the intensity fluctuatio ns in the 
upper photosphere. iFossum fc CarlssonI (|2006f) also used 
TRACE but deduced the upward acoustic flux through 
comparison with hydrodynamic simulations, finding that 
the energy carried by high-frequency waves is insufficient 
to balance the chro mospheric losses. Thi s analysis was 
recently repeated bv lCarlsson et al.l ()2007[ ) using the So- 
lar Optical Telescope aboard Hinode. 

Chromospheric oscillations have generally been mea- 
sured using the Call HandK and infrared triplet lines. 
Imaging measurements using broad filters in H and K 
are mostly dominated by the photospheric contribution 
within the filter bandpass and do not provide a clear indi- 
cation of chromospheric behavior. Traditional slit spec- 
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troscopy gave evidence of significant oscill atory power 
abov e the acou s tic cutoff fr e auencv (e.g. lEvans et aLl 
19631: iLiul [1971 jCra^ fl978l: IMein fc Schmieded 119811 : 
Deubner fc Fleclj 119901). Th e breakthrough simulation 
bv ICarlsson fc SteinI (119971) of th e Call H spectro- 
gram sequence of ^ Lites et al.l ()1993[ ) demonstrated that 
photospheric oscillations produce chromospheric shocks 
through upward propagation, but the role they play in 
chromosphe ric heating and their energetics are still under 
debate fe.g. IKalkofenil2007D . 

In this Letter we analyze a time series of chro- 
mospheric and photospheric velocity images, obtained 
at high spatial resolution with the Interferometric 
Bidimensional Spectrome ter (IBIS; Cayallini 200f|; 
iReardon fc Cavallinill2008t ). We investigate the scaling 
behavior of high-frequency velocity fluctuations in chro- 
mospheric regions with different magnetic topologies, to 
study the basic properties of chromospheric turbulence. 
In Sec. we describe the observations and discuss the 
velocity power spectra, in Sect. [3] we present the analysis 
of statistical properties of velocity increments, and in in 
Sec m our conclusions are summarized. 

2. OBSERVATIONS 

We analyze a dataset obtained with IBIS at the Dunn 
Solar Telescope of the National Solar Observatory on 
2004 June 2 covering an 80" diameter field in the quiet 
Sun at disk center. These data consist of imaging spec- 
tral scans made through the Fel 709.0 nm (height of for- 
mation ^^200 km above rsoo = 1) and Call 854.2 nm 
(height of formation ~1000 km) lines with a spectral 
sampling of 32 and 80 mA, respectively. The overall 
cadence for the complete scan was 19 seconds and the 
scan was repeated 175 times for a total period of 55 min- 
utes. To eliminate a significant source of noise in the as- 
sembled spectra, the narrowband images were carefully 
destretched using broadband images obtained simulta- 
neously with a common shutt er. More detai l is giv en by 
iJanssen fc Cauzzil ()2006[ ) and lVecchio et al.l (|2007[ ). 

Doppler velocities were measured by fitting a second- 
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Fig. 1. — Line-center intensity in the Call 854.2 averaged over 
the duration of the observations. The contours outline the masks 
for the network [blue dotted contours), fibril {green dashed con- 
tours), and internetwork (red solid contours) regions. The substan- 
tial unclassified regions, primarily between fibril and internetwork 
regions, are not analyzed. 



order polynomial to the spectral positions around the 
core of the line after Fourier interpolation onto a finer 
resolution grid and by determining the wavelength shift 
of its minimum. We have performed tests of the accu- 
racy of this approach by adding random noise, consistent 
with the photon statistics of the data, onto an averaged 
reference spectrum and found that, for this spectral sam- 
pling, the precision of the line core position is better than 
1.4 mA or 60 m s~^ for the Fel 709.0 nm line and better 
than 10 ml or 350 m s^^ for the Call 854.2 line. 

We separate several chromospheric regions in the field 
of view on the basis of differences within the maps of 
Call 854.2 line core intensity and the maps of Fourier 
power of the Call line core velocity in the 3- minute 
(4.2 — 6.6 mHz) and 5-minute regimes (3.0 — 3.6 mHz) 
following iVecchio et all (120071). W e add the power con- 
straint because ICauzzi et al.l ( 20081 ) showed that the tra- 
ditional division into network, internetwork, and inter- 
mediate regio ns based on apparent brightness alone (e.g. 
iKriiger et al.l l200lf ) does not effectively separate areas 
with distinctly different chromospheric dynamics, in par- 
ticular the presence or absence of fibrils in the interme- 
diate brightness class. The three areas we thus define — 
namely, network (high- intensity, high 5-minute power), 
fibril (low- intensity, low 3- minute power), and internet- 
work (low-intensity, high 3-minute power) — are shown 
in Fig. m 

For each pixel within these three selected areas we 
performed a one-dimensional Fourier transform of the 
Doppler velocity of Call 854.2 nm line and averaged 
these over that area. For the photospheric Fel 709.0 
nm line the differences in the power spectra among these 
three masks are not significant, and we calculate a sin- 
gle power spectrum for all three areas. Figure [5] shows 
the resulting power spectra, which display characteris- 
tic behavior with a sharp peak around 3.3 mHz for the 
Fel 709.0 line, whereas the Call 854.2 nm line shows a 
broad plateau over 3—5 mHz in the internetwork, demon- 
strating the well-known predominance of 3-minute oscil- 
lations there. 

In this Letter we emphasize the presence of high- 
frequency tails in the average power spectra of each type 
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Fig. 2. — Power spectra of the chromospheric velocities for the 
different types of areas outlined in Fig. [T] and plotted on both log- 
linear (top) and log-log (bottom) scales. The power spectrum for 
the photospheric Fe I 7090 nm line averaged over all three regions 
(purple dot-dashed) is also shown. The horizontal gray lines show 
the noise level for the chromospheric and photospheric power spec- 
tra as estimated in Sec. [2] The linear fits to the high-frequency 
tails are shown in gray on the log-log plots. 



of region. They lie well above our noise estimations 
from the photon statistics (indicated by the horizontal 
lines in Fig. [2]), nearly out to the Nyquist frequency 
of 26 mHz. Such high-frequency tails were already evi- 
dent in tra ditional spectrographic observations by, e.g., 
Evans et al . (1963, Figs. 4 and 5).l0rral]| (Il966l. F ig. 4) 
Woods fc Cram. (|198lL Fig. 1), iDeubner fc Fleckl ()1990l . 
Fig. 2), but our data show them particularly well due to 
a high signal-to-noise and good statistics from the two- 
dimensional spectroscopy. 

The bottom panel of Fig.[2]plots the same power spec- 
tra on log-log scales. The nearly linear trend of the high- 
frequency tails suggests power-law behavior for all three 
chromospheric areas. We fit a line to these tails over the 
range 7—20 niHz (5 — 15 mHz for the network, where the 
power spectra peaks at lower frequencies) and find slopes 
of -2.4,-2.5, and —3.6 for the network, fibril, and inter- 
network regions respectively. A chi-squared goodness-of- 
fit test using these linear fits shows that the tails of the 
power spectra are well represented as a linear trend. A 
1// noise source would also display a power law distri- 
bution, but with a slope of unity. 

The decreasing density with height results in a propor- 
tional increase in the velocity amplitude and would also 
produce an enhancement of the high-frequency power in 
the chromosphere compared to the photosphere. How- 
ever, a process based on a simple density scaling would 
be expected to operate equally at all frequencies and for 
all chromospheric regions (except for a variation in mag- 
nitude due to differing density gradients). Instead, the 
form of the power spectra changes between the photo- 
sphere and the chromosphere, with the appearance of 
the power law trend at high frequencies. Furthermore, 
the slope of this high-frequency tail varies significantly in 
different chromospheric regions, which implies that the 
chromospheric structures play a role in modifying the 
observed velocities beyond a simple scaling of the under- 
lying photospheric motions. We also note that the com- 
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parisons between pow er spectra fro m lin es arising at dif- 
ferent heights made bv lNoved (|1967| ) and lWoods fc CramI 
()1981h clearly show that the power spectra undergo an 
abrupt change near the base of the chromosphere. 

3. TURBULENCE CHARACTERIZATION; VELOCITY 
INCREMENTS 

Although the chromospheric power spectra shown in 
Fig. [2] are not as would be expected from simple scaling 
based on density differences and classical wave propaga- 
tion, they do appear consistent with an energy cascade 
to smaller scales as found in turbulent flows. We there- 
fore focus on the statistical properties of large amplitude 
velocity fluctuations in an approach that differs intrinsi- 
cally from studying solar atmosphere dynamics in terms 
of propagating waves. 

Power spectra can completely characterize turbulent 
fields only if the parameter fluctuations are distributed 
according t o Gaussian probability density functions 
(PDFs; see iFrischI Il995| ) . It is thus important in tur- 
bulence studies to investigate the statistical properties 
across scale separations i to provide information about 
the presence of coherent structures, such as vortices or 
shocks, at the scale I. For a velocity field v{r) the incre- 
ments can be defined as Av^(r) = v(r -I- £) — v(r) and 
usually the longitudinal velocity increments AtJ^n (r) = 
At)£(r) • iji are considered, as in homogeneous and 
isotropic turbulence the energy flux can be expressed 
in terms of the th ird order mo ment of longitudinal ve- 
locity increments (jFrischI 119951 ) . These are also much 
simpler to measure experimentally if the flow speed in 
the probe frame is much larger than the typical veloc- 
ity fluctuations (e.g., solar wind) which implies that the 
timescales can be transforme d into spatial scales accord- 
ing to the Taylor hypothesis (|Frischlll995t ICarbone et al.l 
12004) . This cannot strictly be assumed in our case, but 
nevertheless we analyze line-of-sight velocity increments 
as a function of time separation, as this approach is the 
closest one we can follow with our observational data 
with respect to the standard analysis of longitudinal ve- 
locity increments in turbulence investigations. By denot- 
ing by v(t) the line of sight velocity in a given spatial po- 
sition of the field of view, we consider increments defined 
by AuT-(i) = v(t + t) — where r is the time separa- 
tion. Figure [3] shows PDFs of the velocity increments for 
different time separations, for the three chromospheric 
regions and the photospheric velocities separately. 

The chromospheric PDFs are in all cases non-Gaussian 
and asymmetric, with a prevalence of negative fluctu- 
ations. The asymmetry increases toward small scales, 
with the negative tail becoming more dominant. In com- 
parison, the photospheric PDFs are nearly symmetric, 
but they also change shape with time separation, being 
nearly Gaussian at large scales but developing tails at 
small scales. The change in the shape of the PDFs ob- 
served in both the chromosphere and the photosphere 
indicates a breakdown of self-similarity, which can be 
attributed to intermittency phenomena associated with 
a turbulent cascade. The tails appear to be more de- 
veloped in the chromospheric network than in the other 
regions. 

A measure of the intermittency of velocity increments 
is given by the flatness (i.e. the ratio of the fourth-order 
moment to the square of the second-order moment of the 



increments). The flatness is 3 for a Gaussian distribution. 
The flatness of chromospheric and photospheric velocity 
increments as a function of the time lag r is shown in Fig. 
[4l The flatness of network velocity increments is always 
larger than it is in the other regions, and, most impor- 
tantly, a significant increase of the flatness at small scales 
is found for the network, indicating that intermittency is 
stronger in these regions. 

In order to have a better statistical characterization of 
the small-scale velocity increments in the network. Fig. 
[5] shows the PDFs of the velocity increments conditioned 
by the sign of the velocities, that is, we calculate the 
PDFs of l^Vr{t) = v{t + t)- v{t) for r = 19 s separately 
for the cases in which v{t + r) and v{t) are both positive 
(inward), both negative (outward), or each with opposite 
sign. 

For the latter case a local minimum is present for in- 
crements close to that are obviously less probable when 
the two velocities have opposite signs. The conditioned 
PDF analysis indicates that the tail at negative fluctua- 
tions in the network is mainly due to the contribution of 
structures with inward velocities or opposite velocities. 

4. CONCLUSIONS 

We performed Fourier analysis of spatially and tempo- 
rally resolved measurements in the Call 854.2 nm line. 
We find that high-frequency velocity fluctuations in the 
lower chromosphere follow a clear power-law behavior, 
suggesting turbulence in the chromospheric plasma. We 
studied the PDFs of the velocity increments as a function 
of the time lag in order to characterize this turbulence 
for each of the three defined chromospheric regions (i.e. 
network, fibril, and internetwork) , which showed stronger 
intermittency in the network areas. 

A high-frequency tail is also seen in the photospheric 
velocity up to a frequency of almost 15 mHz. This is 
limited by the higher noise level relative to the observed 
velocities in the photosphere compared to the chromo- 
sphere due to the strong density decrease with height. 
Our observed photospheric power spectra show a smooth 
monotonic decrease with frequency, without an increase 
or structure above the cutoff freq uency as predicted, for 
example, bv lTheurer et al.l (|1997l ) 

Although it is possible that the high-frequency veloci- 
ties seen in the chromosphere are simply due to the up- 
ward propagation of short-period acoustic waves, we find 
that the chromospheric power spectra and the PDFs are 
significantly different from the photospheric power spec- 
trum and PDF, indicating that the velocities in the chro- 
mosphere at these frequencies are not merely the result of 
upward-propagating waves. As Carlsson & Stein (200^ 
point out, only a small percentage of the power present in 
the photosphere at high-frequencies would be expected to 
survive up to higher layers because of radiative damping. 

We suggest instead that the observed chromospheric 
turbulence is generated by the acoustic shocks that are 
present at this height because of the steep vertical den- 
sity gradients and that the non linear shock processes 
produce the cascade of energy to higher frequencies. 
This provides a mechanism for depositing energy into the 
chromosphere that is more evenly distributed both tem- 
porally and spatially than the shock occurrences them- 
selves. Since shocks are observed to be abundant (with 
different characteristics) in both the network and inter- 




Fig. 3. — PDFs of line-of-sight velocity increments calculated from the Call 854.2 nm line for: (a) the network, (b) the fibril, and (c) 
the internetwork regions. The photospheric PDFs, calculated from the Fel 709.0 nm line over all three regions, are shown in panel (d). 
The solid, dashed, and dot-dashed lines refer to time separations t = 19 s, r = 57 s,and t = 760 s respectively. Increments are normalized 
by subtracting their mean and dividing by their standard deviation. 

modeling of this layer of the solar atmosphere. The 
presence of non linearities and the generation of high- 
frequency fluctuations there suggest that traditional os- 
cillation analysis, based on the phase and coherence be- 
tween the photospheric and chromospheric fluctuations, 
may not properly capture the complexities of the dy- 
namic behavior. Extrapolations o f the photosphe r ic en- 
ergy spectrum to h igher layers (jCranmer et al.l 120071 : 
iVerdini fc Velli|[2"007f ) should take into account the mod- 
ification of_the energy spectrum in the chromosphere 
(|Wentzellll977[ ): this modification is due to processes oc- 
curring close to the acoustic cutoff frequency. 

Finally, the presence of significant mechanical energy 
in the chromosphere with a broad spectrum of time 
scale s may play an im portant role in coronal heat- 
ing ([Porter et al .11199^ . The observed high-frequency 
plasma motions are conveniently generated near the sur- 
face where the plasma (3 is close to 1. This allows 
for efficient conversion of these purely acoustic motions 
into a variety of wave modes, through interaction with 
the pervasive magnetic field, and provides many addi- 
tional routes for the continued outward propagation of 
this acoustic energy. If coronal loops are driven at fre- 
quencies comparable to the Alfvcn crossing time (typi- 
cally 10 — 100 sec), the amplitude of Alfvenic perturba- 
tions may increase resonantly and may efficien tly dissi- 
pate the wave ener gy in an impulsive fashion (jHoUwegl 
Il984t iMilano et an il997t iNigro et al.l l2004f ). Acoustic 
shocks are presumably ubiquitous over the solar surface 
and may therefore provide an additional source of energy 
for heating the coronal plasma. 



Fig. 4. — Flatness of line of sight velocity increments as a func- 
tion of the increment time separation r for the different chromo- 
spheric regions as shown in Fig. [T] - the network {dotted line), 
fibril (dashed line), and internetwork [solid line) regions, as well as 
for the photospheric line-of-sight velocity increments {dash-dotted 
line) . 




Fig. 5. — PDFs of chromospheric line-of-sight velocity incre- 
ments calculated in the network for t = 19 s and conditioned by 
the sign of the velocities; the curves refer to the cases in which 
v{t + t) and v{t) are both inward {dashed curve), both outward 
{dash-dotted curve), and of opposite sign {solid curve). Increments 
are normalized by subtracting their mean and dividing by their 
standard deviation. 



network regions of our data (jVecchio et al.|[2008l ). similar 
turbulent cascades are indeed to be expected in both re- 
gions. The resemblance of the network and fibril power 
spectra may result from the direct connection of the mag- 
netic field from the network into the fibrils. 

The presence of non linear turbulent dissipation in 
the chromosphere is an important consideration in the 
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